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ABSTRACT 

We investigate the baryon fraction in dark matter haloes formed in non-radiative gas- 
dynamical simulations of the ACDM cosmogony. By combining a realisation of the Millen- 
nium Simulation (Springel et al.) with a simulation of a smaller volume focussing on dwarf 
haloes, our study spans five decades in halo mass, from 10 10 h^ 1 M Q to 10 15 h~ Y M Q . We 
find that the baryon fraction within the halo virial radius is typically 90% of the cosmic mean, 
with an rms scatter of 6%, independently of redshift and of halo mass down to the smallest 
resolved haloes. Our results show that, contrary to the proposal of Mo et al. (2005), pre- 
virialisation gravitational heating is unable to prevent the collapse of gas within galactic and 
proto-galactic haloes, and confirm the need for non-gravitational feedback in order to reduce 
the efficiency of gas cooling and star formation in dwarf galaxy haloes. Simulations including 
a simple photo-heating model (where a gas temperature floor of Tfl oor = 2 x 10 4 K is imposed 
from z = 11) confirm earlier suggestions that photoheating can only prevent the collapse of 
baryons in systems with virial temperatures T200 < 2.2 Tfl oor « 4.4 x 10 4 K (corresponding 
to a virial mass of M200 ~ 10 10 hr 1 M and a circular velocity of V200 ~ 35 km s _1 ). 
Photoheating may thus help regulate the formation of dwarf spheroidals and other galaxies 
at the extreme faint-end of the luminosity function, but it cannot, on its own, reconcile the 
abundance of sub-L* galaxies with the vast number of dwarf haloes expected in the ACDM 
cosmogony. The lack of evolution or mass dependence seen in the baryon fraction augurs 
well for X-ray cluster studies that assume a universal and non-evolving baryon fraction to 
place constraints on cosmological parameters. 
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1 INTRODUCTION 

In the current paradigm of cosmic evolution, the ACDM model, 
structures grow from an initially smooth density field into a rich 
network of filaments and haloes. Initially, baryons approximately 
follow the collisionless dark matter, but the two components evolve 
differently in non-linear regions after recombination. In proto- 
galactic haloes, for example, the kinetic energy of the collapse 
is thermalised in shocks by the baryons, but rapid radiative cool- 
ing losses allow for further collapse, leading to the formation of 
dense gaseous discs susceptible to swift transformation into stars 
IWhite & Reesll978t) . 

The efficient gas cooling that accompanies the collapse of pro- 
togalactic haloes at high redshift underlies one of the central ap- 
parent conflicts between hierarchical models of structure formation 
and observation. Indeed. IColel 1 1 99 lh and IWhite & Frenkl 1 1991) 
highlighted that, in the absence of additional physics, hierarchi- 



cal models predict that essentially the entire baryonic content of 
the Universe should have cooled at high redshift and, presumably, 
turned into stars by the present epoch. This is in strong contradic- 
tion with observations, which suggest that stars make up less than 
5% o f the baryons in the Universe (for a review, see lBalogh et alJ 
l200ll) . To avoid this 'cooling catastrophe', models of galaxy forma- 
tion frequently invoke various astrophysical mechanisms that coun- 
teract cooling and, in some cases, reheat cold gas. 

One consequence of the cooling catastrophe is the difficulty 
in reconciling the galaxy luminosity function with the ACDM halo 
mass function . This was recognised in the pioneering work of 
White & Rees 1 1978), and it is now widely accepted that a heat- 
ing mechanism (usually referred to as 'feedback') is required to 
explain the reduced star formation efficiency required to match 
the shallow faint end, as well as the sha rp cut-off at the bright 
end, of the galaxy luminosity function (e.g. Kauffmann et al. 1999; 
Somerville & Primacklll999i:ICole et alJl200d : lBenson et alj|2003 : 
Bower et all200dlCroton et all200fj) 
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The main heating mechanisms are commonly believed to be: 
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(i) photo-heating by the energetic photons that reionised the Uni- 
verse at high redshift, (ii) the thermal and kinetic output of evolving 
stars and supernovae (SNe), and (iii) the energy released by accre- 
tion of matter into the supermassive black holes responsible for 
active galactic nuclei (AGN). The importance of each of these pro- 
cesses varies as a function of halo mass. Whereas (i) is thought 
to curtail star formation in extremely low-mass haloes and sub- 
structures, (ii) is assumed to regulate the star formation history 
of dwarf and normal galaxies. The role of (iii) is expected to be 
most relevant for the formation of giant galaxies, typically ellipti- 
cals found predominantly in the dense cores of galaxy clusters (e.g 
iRichstone et alll998l) . 

Of all these, mechanism (ii) is the best studied, but there is 
still no consensus concerning its role in regulating star formation. 
Although it is energetically possible for SNe to suppress star for- 
mation in galactic haloes (Dekel and Silk 1986, White and Frenk 
1991), it has long been recognised that the ef ficiency r equired for 
this process to be viable is uncomfortably high l Benson et al. 2003; 
iKang et all20o3) . Hydrodynamical simulations, for example, show 
that under normal circumstances SNe are far less efficient at ex- 
pelling mass and reheating the IGM than required by galaxy forma- 
tion models in order to match the faint-end of the luminosity func- 
tion iMac Low & Ferrarall999l : IStrickland & Stevensi2000h . Iden- 
tifying an additional or alternative source of feedback is therefore 
highly desirable. 

In a recent paper. Imo et alj l2005t hereafter M05) proposed 
an elegant alternative, namely, that the collapse of large-scale pan- 
cakes and filaments might heat the gas prior to the assembly of 
low-mass haloes. M05 sketch analytic arguments in support of the 
idea that this 'pre-virialisation' gravitational heating might gener- 
ate enough entropy, at low redshift (z < 2), to inhibit substantially 
the accretion of gas into low-mass haloes. If efficient enough, this 
mechanism may offer a simple and attractive resolution to the cool- 
ing crisis, and a possible explanation of the form of the faint end of 
the galaxy luminosity function without the need to invoke feedback 
from non-gravitational sources, such as photoionisation or SNe. 

This is clearly an intriguing proposition, and we present here 
a suite of cosmological gas-dynamical simulations aimed at assess- 
ing the viability of the pre-virialisation heating hypothesis in the 
ACDM cosmogony. Our simulations allow us to probe a vast range 
of halo masses, from 10 10 ft" 1 M to roughly 10 15 ft" 1 M . The 
simulations assume that baryons evolve as a non-radiative fluid; 
this is a conservative assumption when testing the pre-virialisation 
hypothesis, as radiative losses would only serve to facilitate the col- 
lapse of baryons into protogalactic haloes. 

We also use the same suite of simulations to explore the 
baryon fraction at the opposite end of the halo mass function, i.e., 
in galaxy cluster haloes. Adopting a non-radiative gas approach is 
a reasonable simplification here, since the majority of the intra- 
cluster medium (ICM) has a cooling time that exceeds the age of 
the Universe. The most massive galaxy clusters are of particular 
cosmological interest since their baryon fractions are expected to 
trace accurately the cosmic mean. Indeed, the comparison of clus- 
ter baryon fractions (/(,) with the baryon density parameter (Qb) 
implied by Big-Bang nucleosynthesis calculations provides deci- 
sive evidence for a Universe with (dark) matter density well below 
the critical density for closure l White et al. Il993l) . 

In addition, the apparent redshift dependence of cluster baryon 
fractions can be used to constrain the geometry of the Universe 
lSasaki|[l996f) . its deceleration (Pen 1997), and by extension, the 
dark energy equation of state l Allen et all2004l) . These tests exploit 
the redshift dependence of angular diameter distances and rely on 



cluster baryon fractions being roughly universal and non-evolving 
over the redshift range where they can be observed (typically z < 
1). 

Our simulation suite provides the largest sample of haloes 
with non-radiative gas dynamics reported to date. This allows 
us to investigate the mass dependence, evolution and dispersion 
of halo baryon fractions with unprecedented statistical reliability. 
These results can be used to test critically the viability of the pre- 
virialisation heating hypothesis and to examine the stability and 
evolution of cluster baryon fractions. 

In the following section we describe our simulation suite, and 
in Section 3 we outline our analysis methods and main results. We 
discuss them in Section 4 and conclude with a brief summary in 
Section 5. 



2 SIMULATION DETAILS 

To maximise the dynamic range of our halo sample we analyse 
two simulations of different volumes. One (labelled LOW-MASS) 
is a high-resolution 'zoomed-in' simulation of a relatively small 
volume designed to study low-mass haloes. The other (HIGH- 
MAS S)js_a_gas^d^namical realisation of the Millennium Simula- 
tion l Spri nge! etall f2005). a large volume containing many well 
resolved galaxy clusters. We address numerical convergence issues 
by simulating the collapse of a 'pancake' at varying numerical res- 
olution. And, finally, we test explicitly the robustness of our results 
in the LOW-MASS simulation by re-simulating at much higher res- 
olution a single (DWARF) halo with virial mass 10 10 ft" 1 M . The 
numerical parameters and other details of the simulations are listed 
in TableQ 

We evolve our initial conditions i n all cases usin g the pub- 
licly available parallel code GADGET-2 l Sprinsel 2005), with non- 
radiative 1 baryon physics implemented by smoothed particle hy- 
drodynamics (SPH). 

2.1 high-mass simulation 

Our sample of high-mass haloes is drawn from a non-radiative gas- 
dynamical realisation of the Millennium Simulation. It adopts the 
same displacement field and cosmology as the original simulation, 
but has lower mass resolution: 5 x 10 8 gas and dark matter particles 
within a periodic simulation box of side 500 ft" 1 Mpc. The baryon 
density parameter Qb = 0.045 results in particle masses of m gas = 
3.12 x 10 9 ft" 1 M Q ; for the dark matter we adopt f2 dm = 0.25, 
implying a particle mass of md m = 1-42 x 10 10 ft" 1 M©. We 
choose a comoving gravitational softening of 100 ft -1 kpc until 
z = 3, at which point it is fixed in physical units to 25 ft" 1 kpc. 
Further details and analysis of this simulation shall be presented in 
a forthcoming paper (Pearce et al., in prep). 

2.2 LOW-MASS simulation 

In order to analyse a representative region of the Universe whilst 
still resolving low mass haloes, we simulate the evolution of a 
spherical region of radius 7 ft" Mpc, identified at z = in a 
simulation of a box 100 ft" 1 Mpc on a side with the same cosmo- 
logical parameters as the Millennium Simulation. The sphere was 

1 We avoid the common practice of terming non-radiative simulations 'adi- 
abatic', since they include non-adiabatic shocks. 
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1.65 x 10 6 


8.70 x 10 6 
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5.0 


0.040 


0.25 


100 


127 


DWARF 


1.20 x 10 4 


7.80 x 10 4 


3.4 x 10 5 


0.4 


0.040 


0.30 


35.325 


74 


PANCAKE 


1.20 x 10 6 


7.83 x 10 6 


3.5 x 10 5 


10.0 


0.040 


0.30 


N/A 


145 



Table 1. Parameters for the various simulations. Note that each simulation contained an equal number, N p , of both gas and high-resolution dark matter 
particles. Softening lengths, e, are quoted in comoving coordinates hence we mark the HIGH-MASS simulation with an asterisk to denote that its softening was 
switched to physical units at z = 3. The PANCAKE simulation has no boxsize, L, since it was run with vacuum boundary conditions. 



chosen at random from a sample of spherical regions with mean 
density within 10% of the cosmic mean, and devoid of haloes with 
mass M200 > 10 13 Mq, in order to prevent the region from 
being dominated by a single halo. Randomly placed spheres satisfy 
the density selection criterion slightly less than 10% of the time, 
since the volume is dominated by underdense regions. Approx- 
imately one-third of spheres satisfying the density criterion lack 
haloes more massive than 10 13 h,- 1 M . 

Our initial conditions are generated by resampling the re- 
gion with a greater number of particles and adding additional short 
wavelength perturbations, whilst coarse sampling the external mass 
distribution with multi-mass collisionless particles to reproduce the 
large-scale gravitational field. Our resampling algorithm is based 
upon the procedure outlined by Frenk et al. 1 1996), and is described 
in detail by Pow er et al]l2003|) . Gas is added to the high-resolution 
region by splitting each particle into a dark matter particle and a 
gas particle, with mass ratio given by the adopted baryon and dark 
matter density parameters. This implies that each gas particle has a 
corresponding dark matter 'partner' associated with a unique vol- 
ume element at high redshift, a useful feature when tracing the dif- 
ferences in the evolution of the two components in the non-linear 
regime. 

The LOW-MASS simulation features 2.5 x 10 gas particles 
of masses m gas = 1.65 x 10 6 /i -1 M@ and an equal number 
of high resolution dark matter particles, of mass m^m = 8.70 x 
10 6 Mq. At this resolution, the simulated volume yields a 
sample of ~ 1300 well resolved (i.e., A^m > 150) low mass 
haloes at z — whilst remaining relatively computationally in- 
expensive. 

2.3 dwarf simulation 

In order to assess the robustness of our results for low-mass 
haloes we re-simulate a sing le dwarf halo (10 10 /i _1 M ) at much 
higher resolution than its counterparts in LOW-MASS. Because pre- 
virialisation heating is expected to be most effective in haloes as- 
sembling late, we select for resimulation a dwarf halo with rela- 
tively late formation time for its mass. The most massive progenitor 
first exceeds half the final mass of the h alo at z = 0.6, whi lst the 
extensions to the Press-Schechter t heory IPress & SchechteJl974t) 
described bv lLacev & Cold 1 19931) suggest that the most probable 
formation time for a halo of this mass is z ~ 2. 

We selected the halo from a parent simulation of box length 
35.325 Mpc and density parameters (fi m ,fiA = 0.3,0.7). 
We apply the same resimulation technique as for LOW-MASS, us- 
ing 3.4 x 10 5 gas and high-resolution dark matter particles. We 
adopt a baryon density parameter of 0,^ = 0.04, which implies 
particle masses of m gas = 1.20 x 10 4 h" 1 Mq and md m = 
7.8 x 10 4 ft -1 M Q . At z = 0, the halo has mass M 20 o = 



9.5 x 10 h Mq, virial radius r2oo = 34.5 h kpc, and cir- 
cular velocity V200 = 35 km s -1 . 



2.4 pancake simulation 

The pre-virialisation mechanism outlined by M05 assumes that gas 
is heated during the pancake-like collapse of the large-scale struc- 
ture within which dwarf haloes are embedded. It is therefore im- 
portant to explore whether our numerical techniques are suitable to 
describe this process accurately, as well as whether the results are 
not artificially marred by limited numerical resolution. 

We investigate this by simulating the collapse of an idealised 
pancake with similar numerical resolution as that of the LOW-MASS 
simulation, and then vary the resolution in order to assess conver- 
gence. The simulation involves the collapse of a uniform spherical 
region of mass 3 x 10 12 ft -1 Mq, initially perturbed so that it col- 
lapses along one axis to form a flattened pancake at z = 2. We 
use a 'glass' rather than a grid in order to minimise the artifacts 
introduced by anisotropies in the grid. The desired dynamics are 
achieved by compressing the sphere along one axis and expanding 
it by half of the compression factor along the other two axes, with 
initial velocities computed using linear theory. 

The choice of parameters for this simulation is motivated by 
the discussion of M05, who argue that such 'pancakes' might repre- 
sent the typical environment where dwarf galaxy halos are formed, 
and that pre-virialisation heating might have been missed in early 
simulations because of inadequate resolution. As noted by M05, 
a collapsing pancake forms shocks on both sides, so a faithful 
treatment of its thermodynamic evolution requires > 8 smoothing 
lengths across its collapsing axis. Following M05, and assuming 
that shocks operate as the collapsed structure approaches an over- 
density of ~ 10, this implies a pancake thickness of 200 h~ kpc 
(comoving) and a minimum resolution of 25 h~ kpc, again co- 
moving. We fulfil this criterion using a particle mass similar to that 
of the LOW-MASS simulation. At an overdensity of 10 QbPc(z) a 
mass of 40 m p particles (comparable to what is used to define the 
SPH smoothing length scale) is contained within a sphere of co- 
moving smoothing length h Brn i ~ 25 h~ kpc. 



3 ANALYSIS & RESULTS 

Figure Q shows various snapshots of the LOW-MASS simulation. 
Each box is 10 h^ 1 Mpc on a side and shows the positions of 
gas particles, colour-coded according to density. This figure shows 
clearly the highly anisotropic nature of the large-scale structure, but 
also highlights the points that many dwarf dark matter halos had al- 
ready collapsed by z ~ 2. As we discuss below, this has important 
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Figure 1. Redshift-progression of the projected gas density within a cube of side 10 h^ 1 Mpc (comoving) in the LOW-MASS simulation. The progression 
clearly illustrates that by z = 2 many dark matter haloes have already collapsed, driving their associated gas to high overdensities. Such gas is afforded 
stability against the shocks that develop as the large-scale environment, within which it is embedded, collapses into pancakes and filaments 



implications for the efficiency of pre-virialisation heating and its lated region at z 
effects on the baryon fraction of collapsed systems. effects. 



0, since they may have been subject to boundary 



3.1 Halo finding algorithm 

We u se a friends-of-friends group finding algorithm jDavis et"all 
1985) with a short linking length (6 = 0.05) to locate the cores of 
haloes in the simulation volumes. The centres of these haloes are 
used as starting points for an iterative algorithm that finds spher- 
ically overdense regions of mean enclosed density 200p C rit (z), 
whose centres coincide with their centres of mass. We define the 
virial radius, r2oo, as the radius of this sphere. Other 'virial' quan- 
tities quoted for haloes refer to measurements within this radius, 
unless otherwise specified. 

We consider only haloes with at least 150 dark matter parti- 
cles in order to minimise the effects of poor numerical resolution, 
and clean the sample by removing haloes partially contained within 
other haloes. In the case of LOW-MASS we also disregard haloes 
located closer than 200 h" 1 kpc from the boundary of the resimu- 



3.2 Baryon fractions 

3.2.1 HIGH-MASS 

The right-hand panels of Figure [2] show the baryon fraction in 
haloes identified in the HIGH-MASS simulation at z = 1 and at 
z = 0. These panels show the largest sample of haloes simulated 
with non-radiative gas physics reported to date, with approximately 
49, 500 and 115, 000 haloes resolved at z = 1 and 2 = respec- 
tively. The panels show that the baryon fraction is independent of 
mass, and has not evolved at least since z = 1, the highest redshift 
for which cluster baryon fractions can be estimated reliably from 
X-ray observations (e.g. the observations of Allen et al. ( 2004j) span 
the redshift range 0.07 < z < 0.9). 

The mean cluster baryon fraction within the virial radius is 
approximately 90% of the cosmic mean, with relatively small scat- 
ter; the root-mean-square dispersion is less than 3% for haloes of 



The effects of pre-virialization and photoheating on the baryon fraction of ACDM haloes 5 



'°9lO N drr. 

2.0 3.0 4.0 5.0 2.5 3.0 3.5 4.0 4,5 5.0 




9.Q 1Q.Q 11.0 12.0 12.& 13.0 13.5 14.0 14.fi 1 5.D 15.5 



lo 9io W 2oo [ h 1|J s] 

Figure 2. Baryon fractions, in units of the universal value, of well resolved haloes (N^m > 150) drawn from the LOW-MASS (left-hand panels) and HIGH- 
MASS (right-hand panels) simulations at z = 1 (upper panels) and z = (lower panels). The large dots and error bars show the mean and rms of the 
distribution, respectively. The green triangle marks the baryon fraction of the halo in the DWARF simulation. The upper horizontal axis gives the equivalent 
number of dark matter particles at a given mass scale. The downward arrows in each plot illustrate the mass scale corresponding to 500 dark matter particles. 



mass M200 > 3 x 10 14 hT 1 M©, and the difference between the 
10 th and 90 th percentiles is always less than 7.5% over the same 
mass range. The scatter remains small for all well resolved haloes; 
the rms scatter is less than 6% for haloes resolved by at least 500 
dark matter particles. This result is in broad agreement with pre- 
vious simulations of cluster baryon fractions in the non-ra diative 
regime using SP H codes (e.g.lNavarro et alll995tlEke et alll998t 
iFrenk et alll999tlEttori et alfcOOd) . ~ 

These results are only weakly dependent on radius; within a 
radius encompassing a mean inner density 500 times greater than 
critical, rsoo, the results are very similar to those within V2oo' for 
M200 > 3 x 10 14 ft -1 M© the baryon fraction remains, on av- 
erage, 90% with an rms dispersi on of 3%. This radius, w hich is 
~ 0.7r 20 o for an NFW profile iNavarro et ai]ll99d 1 19971) with 
concentration c = 5, is roughly the maximum radius for which 
total cluster masses can be estim ated reliably from X-ray observa- 
tions fe.g. lVikhlinin et all20o"r3) . 



3.2.2 LOW-MASS 

The left-hand panels of Figure [2] show that a similar result applies 
to the halo sample identified in the LOW-MASS simulation. The 
baryon fraction of galactic and dwarf haloes is also about 90% of 
the cosmic mean, and shows no discernible dependence on mass or 
redshift. 

A noticeable downturn is observed below about 10 10 h~ 1 M@, 
but this may be ascribed to the poorer resolution affecting such 
haloes, since an underestimate of the gas density at accretion 



shocks leads to artificially high post-shock entropies. Note that a 
similar downturn is also seen in the HIGH-MASS sample for haloes 
resolved with fewer than ~ 500 particles. Downward arrows show 
the mass scale corresponding to 500 dark matter particles in each 
simulation. 

3.2.3 DWARF 

Figure|2|suggests that for haloes resolved with more than ~ 500 — 
1000 particles, numerical resolution does not affect the measured 
baryon fractions. Further supporting evidence comes from the re- 
sults obtained for the high resolution individual DWARF halo. These 
are shown in Figure|2|with a filled triangle, and are consistent with 
LOW-MASS haloes of similar mass. 

3.2.4 PANCAKE 

One may still worry that, despite the apparent convergence of the 
results shown in Figure[2] the numerical resolution is insufficient to 
capture the shocks during the pancake collapse phase that accom- 
panies the formation of haloes. We can test this by examining the 
PANCAKE simulation. 

Figure[3]shows the evolution of the entropy (measured by s = 
T/nV 3 ), and of the temperature of the gas during the collapse of 
the sphere to a pancake configuration. Note that when computing 
the electron density, n e , we assume the gas is fully ionised and of 
primordial composition. At z — 2 the system reaches maximum 
asphericity, and is well described by a plane of comoving thickness 




Figure 3. The evolution of the median entropy parameter ( top) and median 
temperature (bottom) of gas particles within the PANCAKE collapse simula- 
tion. The solid line in each case represents the simulation featuring particle 
mass resolution identical to the LOW-MASS simulation. We supplement this 
run with degraded resolution runs with particle number reduced by a factor 
of 8 (dotted line) and 64 (dashed line). An improved resolution run is also 
presented, with particle number a factor of 8 greater (dot-dashed line). The 
vertical dotted line marks z = 2. 



Figure 4. Baryon fractions, in units of the universal value, at z = for 
the LOW-MASS simulation with and without our photoionising background 
model. The open and solid green triangles show the baryon fraction of the 
DWARF halo with and without the same model, whilst the vertical dot- 
dashed line marks the virial mass and virial temperature equating to the 
temperature floor. The solid line denotes the baryon fractions predicted by 
M05 for haloes in the presence of a UV background, whilst the dotted line 
shows the prediction of the M05 gravitational preheating model. 



25 h' 1 kpc (c.f. comoving softening of 10 h 1 kpc) and comoving 
radial extension of 1.5 /i -1 Mpc. 

The collapse heats gas to a median temperature of 3 x 10 4 K, 
and the median specific entropy reaches s ~ 1 keV cm 2 . We test 
explicitly the effect of resolution by re-running the PANCAKE simu- 
lation with particle numbers reduced by factors of 8 and 64 relative 
to the standard resolution of LOW-MASS; we complete the resolu- 
tion study by running a higher-resolution case where the number 
of particles was increased by a factor of 8. Gravitational softenings 
were scaled as e oc Np^ 3 . 

As shown in Figure|5]the post-shock median entropy and tem- 
perature jumps are quite insensitive to numerical resolution, al- 
though the transition becomes noticeably sharper as the resolution 
increases. This test shows that limited resolution does not lead to a 
substantial underestimate of the entropy jump in numerical simula- 
tions of pancake-like collapse. Actually, poor resolution leads typ- 
ically to underestimation of the true densities: entropies are there- 
fore typically overestimated in poor resolution simulations. This 
result, together with the consistency between the LOW-MASS and 
DWARF simulations, gives us confidence that baryon fractions in 
our simulations are not unduly affected by resolution effects. 

3.3 Photoheating 

The one feedback mechanism that is certainly present at early times 
is associated with the energetic photons that reionised the Universe 
at high redshift. This has long been recognised as having the po- 
tential to inhibit the formation of galaxies in low-mass haloes, al- 
though there is still no consensus concerning the mass scale be- 
low which photoheating becomes effective at halting galaxy forma- 
tion jBlumenthal et alJll984t lEfstathioulll992t buinn et alj|l996t 
iBullock et all200 0: Benson et al. 2002f)7Most semianalyti c galaxy 
formation models have so far adopted the prescription of Gnedin 
( 2000) to determine the gas accreted by haloes for given IGM pres- 



sure, but recent results presented by H oeftetalj f2006) suggest that 
Gnedin's approach may substantially overestimate the mass scale 
of photoheating. 

This motivates us to include a simple photionisation heating 
model in our simulations. The aim is twofold. On the one hand, we 
wish to shed light on the disagreement about the effects of photo- 
heating on gas fractions, but, on the other hand, we would also like 
to explore whether photoheating may act to suppress the efficiency 
of gas accretion in the early phases of the hierarchy, facilitating and 
enhancing the thermodynamic effect of pancake-driven shocks. 

We investigate the combined effect of gravitational and pho- 
toheating by re-running the DWARF and LOW-MASS simulations, 
again with non-ra diative gas physics, but, motivated by the third- 
year WMAP data iSpergel et ai]|200d) . imposing a spatially uni- 
form temperature floor for all gas particles at z = 11. To aid 
as much as possible the pre-virialisation generation of entropy, 
we adopt for the temperature floor a rather high value, Tfl oor = 
2 x 10 4 K, consistent with the maximum temperature of the IGM at 
mean density, as probed by QSO absorption spectra I Scha ve et alj 
2000). 

This should clearly impact gas accretion on dwarf galaxy 
haloes and, in particular, our DWARF halo, where the virial tem- 
perature is only ~ 4.4 x 10 4 K at z = 0, only a factor of 2.2 above 
Tfloor- We therefore expect that a considerable fraction of the gas 
bound to small halo progenitors at high redshift should be photo- 
evaporated from this structure. 

Figure|4]illustrates the effect of the additional heating on the 
baryon fractions of the LOW-MASS halo sample (open circles) and 
compares them with the results of the non-radiative run (filled 
circles). Photoheating introduces a well-defined mass scale be- 
low which gas accretion is strongly suppressed. Below A/200 ~ 
10 10 M© haloes are able to retain less than one half of their 
share of baryons within their virial radii; the effect is as large as 
90% in haloes below 3 x 10 9 ft -1 Mq, corresponding to an effec- 
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Figure 5. A comparison of the evolution of the mass of the main DWARF 
halo progenitor (black solid line), with the evolution of the final halo gas 
mass that is already collapsed at redshift z, as defined by the condition that 
its density exceeds 10p c (z) (green dotted line). The data are normalised to 
their values at z = 0. The latter quantity is also shown for the case when our 
photoheating treatment is applied (red dashed line), this time normalised to 
the 2 = value of the purely non-radiative case. 



tive virial temperature very similar to the photoheating temperature 
floor. 



4 DISCUSSION 

The main result of the previous section is that, in the non-radiative 
approximation, the baryon fraction of ACDM haloes is independent 
of redshift as well as of mass, in the range resolved by our simu- 
lations (10 10 -10 15 /i -1 Mq), Photoionisation reduces the baryon 
fraction only in haloes with virial temperature comparable to that 
imposed by the ionising photons, typically just above 10 4 K. We 
discuss below the implication of these results for models of galaxy 
formation and for the interpretation of observations of baryon frac- 
tions in clusters. 



4.1 Baryon fraction bias 

The lack of dependence of baryon fractions on halo mass is in- 
triguing, as is the fact that the mean value within the virial radius 
is only ~ 90% of the cosmi c mean. The same result has been ob- 
served in other simulations iNavarro et alll 19951 lEke et alj[l99l 
Fren k et alj 1999), and has been ascribed to the collisional vs. col- 
lisionless nature of the baryons and dark matter, coupled to the hi- 
erarchical assembly of haloes in the ACDM cosmogony. Indeed, 
during the many mergers that mark the formation of a halo, shocks 
act to stop the gas whilst the dark matter streams through freely. 
This leads to a temporary spatial offset between dark and gaseous 
components during which energy and angular momentum are trans- 
ferred from the dark matter to the baryons, as discussed in detail by 
lNavaffo&WhItel<1993l) . The energy gained during mergers results 
in a more extended gaseous component, and an overall (slight) re- 
duction in the baryon fraction relative to the cosmic mean. 



4.2 Pre-virialisation heating 

The results of Figure[2]imply that pre-virialisation heating is inef- 
fective at preventing the collapse of baryons into low-mass haloes, 
even for masses as low as 10 10 /i -1 Mq. Figure [4] compares the 
baryon fractions of LOW-MASS haloes at z = (filled circles) with 
the predictions of the M05 model (dotted line). M05 argue that 
the IGM in low-mass haloes should have been heated by shocks to 
roughly ~ 10 kev cm 2 prior to halo assembly, and that this would 
lead to a reduction of > 50% in the baryons filling haloes of mass 
< 6 x 10 11 /i" 1 Mq . M05 propose a fitting formula to characterise 
this effect; 

/■halo -i 

— = " CD 

ycosmic (l + M c /M) a 

where a = 1 and M c = 5 x 10 11 h~ x Mq is a characteristic 
mass scale. This function clearly fails to reproduce our results, and 
suggests that the hypotheses on which M05 base their model are 
not satisfied in our simulations. 

The main premise of M05's model is that most low-mass 
haloes form in extremely aspherical regions where their assembly 
might be delayed, allowing for pancake-driven shocks to elevate 
the entropy of the IGM prior to halo assembly. Our simulations, 
however, indicate otherwise. 

Firstly, low-mass haloes surviving to the present were, at the 
time of their formation, in regions where pre-virialisation shocks 
were weaker than envisaged by M05. For example, we find a typi- 
cal post-collapse entropy of ~ 1 keV cm 2 in our pancake collapse 
simulations, about an order of magnitude lower than adopted by 
M05 to compute the model shown by the dotted curve in Figure|4] 
These results suggest that the halo assembly process is approxi- 
mately scale-free; if pre-virialisation heating does indeed occur, it 
affects all haloes in similar measure, leaving no particular signature 
in low-mass haloes. 

Secondly, the material destined to make a low-mass halo col- 
lects into dense, early-collapsing clumps prior to the collapse of 
the surrounding 'pancake'. The pancake, in other words, is not a 
nearly uniform aspherical structure where shocks may propagate 
freely, but rather a large-scale feature where a substantial fraction 
of the mass is in collapsed clumps. 

We illustrate this in Figure|5| where we plot (dotted line) the 
fraction of the final DWARF halo gas that resides in collapsed struc- 
tures, as quantified by the condition p > 10 p c tit(z). By this rather 
strict measure, half the DWARF gas is already in collapsed structures 
prior to z ~ 2, although by then the most massive halo progenitor 
(solid line) has only about ~ 25% of the final mass. This early 
aggregation of the halo gas into dense structures prevents it from 
being shock-heated by the pancake-driven shocks, reducing further 
the pre-virialisation heating efficiency. 

4.3 Photoheating 

As shown in Figure|4] the baryon fraction may be reduced because 
of heating by a photoionising background, but the effect (at z = 0) 
is restricted to haloes with virial temperatures < 2.2Tfl oor . This 
implies that ionising photons, which are unlikely to heat the gas 
to temperatures much higher than ~ 2 x 10 4 K are only able to 
influence the formation of galaxies in haloes less massive than ~ 
10 10 /i -1 Mq. Our results for the baryon fraction in this case are 
well described by Eqn.Q but with the revised parameters suggested 
by M05 for their photoheating model: a — 3 and M c — 1.7 x 
10 9 ft -1 M©; we show this fit in Figure[4]as a solid line. 
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The results shown in Figure|4]agree with those of lHoeft et alj 

( 2006), who used a more detailed treatment of the UV background; 
the similarity of our findings is rather encouraging. We concur with 
their assessment that the characteristic mass scale of photoevapo- 
ration is probably considerably lower than derived from the filter- 
ing mass formalism o flGnedml feoOO). Thus, although photoheating 
can reduce the baryon fraction in low-mass systems, it appears to 
be less efficient at shaping the extreme faint end of the galaxy lu- 
minosity function than previously inferred through semi-analytic 
modelling (e.g. lBenson et all2002tlsomervillel2002h . 

Our results seem to be robust to numerical resolution, as 
shown by the good agreement between the baryon fraction of the 
LOW-MASS and the DWARF haloes plotted in Figure^ In the latter 
case, photoheating reduces the baryon fraction by 50% but there is 
no evidence that pre-virialisation has played a role; indeed, inspec- 
tion of the evolution of the collapsed gas fraction (red dashed line 
in Figure |5J shows no discernible feature that may be associated 
with the collapse of the large-scale structure. 

M05's pre-virialization model, shown as a dotted line in Fig- 
ure|4] requires baryon fractions to be halved in haloes as massive as 
5 x 10 11 M Q ;this is an order of magnitude larger than the mass of 
haloes significantly affected by photoheating and pre-virialisation 
in our simulations. Our results suggest, then, that other feedback 
mechanisms are required to match the faint end of the galaxy lumi- 
nosity function in the ACDM scenario. 



4.4 Application to cluster surveys 

Our HIGH-MASS simulation also features over 115,000 galaxy 
cluster-sized haloes, and demonstrates that, in the non-radiative 
regime, cluster baryon fractions are independent of virial mass, dis- 
play little dispersion, and do not evolve significantly over the red- 
shift range < z < 1. Observationally, mass profiles of clusters 
are typically only estimated reliably out to a maximum radius of 
rsoo; the results within that radius are very similar to those plotted 
in Figure|2|at the virial radius. In haloes where the region interior 
to rsoo is resolved by at least 500 dark matter particles, the mean 
cluster baryon fraction within rsoo at z = remains approximately 
90% of the cosmic mean, again with an rms dispersion of ~ 6%. 

The key to the applicability of our results to cosmological tests 
is the validity of our non-radiative treatment of the ICM. Whilst the 
failure of purely non-radiative models to match some of the global 
scalin g relations exhibited by clusters is w ell documente d (e.g. 
[ Evrard & HenrJl99ll;rKaiserll99ll;lNavarro et alll995l ; lEke et alJ 
1998), this does not necessarily imply that non-radiative models 
give the wrong fraction of hot baryons in observed clusters. Since 
the p 2 dependence of thermal bremsstrahlung implies that the X- 
ray emissivity of clusters is dominated by the central region, it is 
possible to obtain agreement with the observed X-ray scaling rela- 
tions by modifying only the central gas density (for instance, with 
radiative cooling or AGN feedback), whilst leaving the density of 
the bulk of the gas un changed (e.g. Balog h et alJll999tlVoit et alJ 
l2002ll2003tlMcCarthv et all200ll20Ml 

It should be noted, however, that the X-ray luminosity- 
temperature relation can also be explanied by models that do af- 
fect baryons at large radii (e.g. iKav et af] 12004 iKravtsov et alJ 
l2005t Ettori et al. 2006). In such models the hot baryon fraction 
at rsoo can be reduced by up to ~ 30%, even for clusters with 
M200 > 10 15 h~ x Mq. However, such large-scale reductions 
by means of cooling and star formation conflict with optical con- 
straints (e.g. Balogh et alj|200ll) . and the necessary level of feed- 



back from SNe also app ears unfeasibly high ( Benso n~et alJl2003t 
IScannapieco & Ohl2004 . 

The effect of non-gravitational processes on the large-scale 
properties of rich clusters remains a source of debate. We antic- 
ipate that analyses of the large samples of rich clusters provided 
by the Chandra and XMM-Newton observatories will foster the de- 
velopment of a definitive picture of the ICM. This picture may well 
show that baryons at large radii are only minimally affected by non- 
gravitational processes, thus validating our results that rich cluster 
baryon fractions do not evolve for z < 1 and exhibit little disper- 
sion. 



5 CONCLUSIONS 

We have measured the baryon fractions of a large sample of haloes 
drawn from a suite of non-radiative gas-dynamical simulations of 
the ACDM cosmology. The haloes span five orders of magnitude in 
virial mass, from dwarf galaxy haloes to large clusters. Within the 
virial radius, the baryon fraction averages 90% of the cosmic mean, 
with a fairly small dispersion (~ 6% rms) and shows no depen- 
dence on redshift for well-resolved systems. This is at odds with 
the 'p re-virialisation' gravitational heating proposed by M o et alj 
1 2005). Pre-virialisation, if at all present, plays only a minor role in 
setting the budget of baryons that accrete into low-mass haloes. 

Photoheating, modelled here as resulting from a uniform tem- 
perature 'floor' of 2 x 10 4 K imposed on the baryons from z = 11, 
is only able to reduce the baryon fraction in haloes with virial tem- 
peratures comparable to the photoheating floor. The absence of a 
strong mass trend in the baryon fractions of low mass haloes high- 
lights the need for non-gravitational feedback as a means to regu- 
late gas cooling and star formation in low-mass haloes, in order to 
reconcile the ACDM halo mass function with the observed galaxy 
luminosity function. 

At rsoo, the typical maximum radius at which current X-ray 
observatories can probe cluster temperatures, the baryon fraction 
remains similar to that at r2oo, again with similarly small disper- 
sion. Provided the effects of non-gravitational physics within rsoo 
on rich cluster baryon fractions are small, we conclude that mea- 
surements of the matter density parameter and the dark energy 
equation of state based on cluster baryon fractions should not suf- 
fer from large uncertainties arising from the dispersion of baryon 
fractions, nor from redshift dependent systematics. 
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